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Abstract 



We present calculations of next-to-leading order QCD splitting functions, employing the 
light-cone gauge method of Curci, Furmanski, and Petronzio (CFP). In contrast to the 
'principal-value' prescription used in the original CFP paper for dealing with the poles of 
the hght-cone gauge gluon propagator, we adopt the Mandelstam-Leibbrandt prescription 
which is known to have a solid field-theoretical foundation. We find that indeed the 
calculation using this prescription is conceptionally clear and avoids the somewhat dubious 
manipulations of the spurious poles required when the principal-value method is applied. 
We reproduce the well-known results for the flavour non-singlet splitting function and the 
N^, part of the gluon-to-gluon singlet splitting function, which are the most complicated 
ones, and which provide an exhaustive test of the ML prescription. We also discuss in 
some detail the x = 1 endpoint contributions to the splitting functions. 



1 Introduction 



The advantages of working in axial gauges when performing perturbative QCD calcula- 
tions are known since a long time Those gauges enable us to retain, in higher order 
calculations, a natural 'partonic' interpretation for the vector field, typical to leading log 
approximation. 

Among axial gauges, the one which enjoys a privileged status is the light-cone axial 
gauge (LCA), characterized by the condition n^Af^ = 0, being a light-like vector 
(n^ = 0). At variance with temporal {n^ > 0) and spacelike (n^ < 0) axial gauges, 
which do have problems already at the free level 0, and with the spacelike 'planar' 
gauge 0] in which the behaviour of the theory in higher loop orders is still unsettled 0], 
LCA can be canonically quantized |Q and renormalized ||^ at all orders in the loop 
expansion following a well-established procedure. To reach this goal it is crucial to treat 
the 'spurious' singularity occurring in the tensorial part of the vector propagator 

according to a prescription independently suggested by Mandelstam p and Leibbrandt 
(ML) and derived in Ref. in the context of equal-time canonical quantization: 

J-^J-^ I = "•' . (2) 

(nl) [nl] nl + iesign(n*/) n*lnl + ie 

the two expressions being equal in the sense of the theory of distributions. The vector n* 
is light-like and such that n*n = 1. Denoting by /_l the transverse part of the vector 
orthogonal to both and n* , one has 

2{nl){n*l) =1^ + ll . (3) 



The key feature of the ML prescription is that the spurious poles in the complex Iq 
plane are placed in the second and fourth quadrants, i.e., with the same pattern as one 
encounters for usual covariant denominators like + ie). One can therefore perform 
a proper Wick rotation to Euclidean momenta, and a suitable power counting criterion 
in the Euclidean integrals will give information on the ultraviolet (UV) divergencies of 
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the corresponding Minkowskian integrals. This is in contrast to the Cauchy principal 
value (PV) prescription, which, under a Wick rotation, entails further contributions and 
therefrom a violation of power counting. 



A crucial property of the ML distribution is the occurrence of two contributions with 
opposite signs in the absorptive part of the vector propag ator [0, 

disc[P,.(/)] = 2vr^(/^)e(/o) f-^7.. + ^- ''"^%t'''^^" ) (4) 

-27r + Z^)e(/o)^^^^^^^^^tZ!fll^i . 

n*n f]_ 

Here the second, ghost-like, contribution (which is not present in the PV prescription) is 
responsible for the milder infrared (IR) behaviour of the ML propagator. The presence of 
this axial ghost was stressed in [||; its properties are exhaustively discussed in 0. Clearly, 
if one has a cut diagram with, say, m final-state gluons, there is a discontinuity like (^) 
for each of the gluons, i.e., the phase space will split up into 2™ pieces. 

One of the most interesting and non-trivial applications of the LCA is the computa- 
tion of the (spin independent) splitting functions for the two-loop Altarelli-Parisi (AP) 
evolution of parton densities, following a method proposed and used by Curci, Furmanski 
and Petronzio (CFP) in Refs. TD|. This method is based on the observation [rT| that in 



axial gauges the two-particle irreducible kernels of the ladder diagrams are finite, so that 
the coUinear singularities that give rise to parton evolution, only originate from the lines 
connecting the kernels. Therefrom, using renormalization group techniques, the splitting 
functions are obtained by some suitable projection of the ladder diagrams, exploiting the 



factorization theorem of mass singularities We refrain from giving further details of 



the CFP method, since these can be found in We just mention at this point that 

one projects on the quantity Fjj, given by 



Fjj I X, O-si ^ j 



5(1 - x)bi. +xFF [ -^5 (x-—] UiK L, 

_ ' ^ J {2-kY V pnj 1-VK ^ 



(5) 



where 2e = (4 — d) and is a 2PI kernel, which is finite in the light-cone gauge 0. 
The labels i,j run over quarks and gluons; in the flavour non-singlet case one has just 
{ij} = {qq}. Furthermore, in (||) PP denotes the pole part, and the projectors Ui,Li are 
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given by 



d~2 



^ pn 



(6) 



The splitting functions Pij to the desired order can be read off from the 1/e pole of Fy-: 

r,,(.,„.,i)=i(i-.)i,-l(gp,fW + i(|)'n™W + ...)+o(i) . (7) 



For future reference, we write down a similar expression for the residue Zj {j = q, g) 
of the pole of the full quark (or gluon) propagator: 



Inspecting Eqs. , (0) , (§|) , we see that Zg and Zg contribute to the endpoint (~ 5(1 — x)) 
parts of the splitting functions Pgg and Pgg, respectively. 



In the above references [y, |T0[, the 'spurious' singularity 1/nl of the gluon propa- 
gator was handled according to the PV prescription. The method of |iy] has been 
very successful in providing the first correct result for the next-to-leading order (NLO) 
gluon-to-gluon splitting function. The result previously obtained in the operator product 



expansion (OPE) method was not correct due to a subtle conceptual problem which 
was recently clarified |16[. The new Feynman gauge OPE calculations confirmed 
the old CFP result. Despite of this success, the LCA calculation with PV prescription 
is considered dubious because of the difficulties with power counting and Wick rotation 
mentioned above. In particular it is not clear whether its 'calculational rules' remain valid 
in higher orders. We note that the precise description of some of the new high precision 
collider data call for the extension of the NLO QCD analysis to next-to-next-to-leading 
order (NNLO). Therefore a deeper understanding of the formal field-theoretical basis of 
the CFP method is strongly motivated. The use of PV is by no means mandatory in 
the CFP method, it can also be applied when handling the 1/ln singularities with the 
theoretically more sound ML prescription. 



A first attempt using the ML prescription in connection with the CFP method has 
been performed in Ref. |^, where the one- loop AP splitting functions [|T^ have been 
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correctly reproduced, both for the flavour non-singlet and for the flavour singlet case. A 
new characteristic feature of this calculation is that 'real' and 'virtual' contributions are 
separately well-defined in the limit x ^ 1, x being the longitudinal momentum fraction, 
at variance with the corresponding PV result. This occurs thanks to the presence of the 
'axial' ghost, which, standing by the usual gluon term, protects its singular behaviour 
with respect to the transverse momentum. There is no need of any IR cutoff to regularize 
intermediate results. 

Beyond one loop, the calculation of the splitting functions according to the CFP 



method in LCA with the ML prescription, has already been tackled in a recent paper |T7 . 
We believe, however, that improvements to the calculation can and should be made. 
First of all, only the Cp part of the flavour non-singlet splitting function is studied 
in [|l^. In this paper we will also calculate the CpTf part and, in particular, the far more 
complicated piece ~ CpNc of this function, as well as the Nq part of the gluon-to-gluon 
splitting function contributing to the flavour singlet sector. As we will show, this set of 
functions we consider comprises all possible one-loop structures of QCD and thus enables 
an exhaustive test of the ML prescription in this application. The ML calculation of the 
other singlet splitting functions, like the non-diagonal quark-to-gluon (and vice versa) 
one, is therefore not really required in this context: they will certainly come out correctly 
if the prescription works for the far more complicated cases we study. 



Secondly, the power and virtues of the ML prescription were not fully exploited in , 
where some contributions resulting from the axial ghosts of the ML prescription were ne- 
glected. These contributions are ~ 6{1 — x) and thus only affect the endpoints of the 
diagonal splitting functions. Nevertheless, their inclusion is required for a complete anal- 
ysis, since only then the crucial question of the finiteness of the two-particle-irreducible 
(2PI) kernels in the light-cone gauge can be fully answered. We also remind the reader in 
this context that in the original CFP papers |TU|, the endpoint contributions to the di- 
agonal splitting functions were never determined by explicit calculation, but were derived 
in an indirect way from fermion number and energy-momentum conservation. The fact 
that we pay more attention to the point x = 1 will enable us to improve this situation 
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to a certain extent: for the first time within the CFP method, we will determine the full 
part ~ CpTf6{l — x) of P^q^^^ by explicit calculation. 

Finally, in |T7] a principal value regularization was still used at some intermediate 



steps of the calculation. Even though this was only done at places where it seemed a safe 
and well-defined procedure, it is more in the spirit of the ML prescription to abandon the 
PV completely and to stick to one single regularization, the dimensional one. This view is 
corroborated by the observation that the PV regularization as used in [|17| actually turns 
out to become technically too complicated when one studies the CpNc part of the flavour 
non-singlet splitting function, or the N^, part of Pgg^ ■ 

The remainder of this paper is organized as follows: to set the framework, we will 
present a brief rederivation of the leading order (LO) quark-to-quark splitting function 
Pqq^ in sec. 2. Section 3 will contain the calculation of the flavour non-singlet splitting 
function at two loops. More specifically, we will discuss in detail the treatment of the 
various virtual-cut and real-cut contributions in subsections 3.1 and 3.2, respectively, 
while sec. 3.3 presents the final results of the calculation. In 3.4, we discuss the endpoint 
contributions and provide a sample calculation of a two-loop contribution to the quark 
self-energy in the LCA with ML prescription. Section 4 deals with the calculation of the 
A^^ part ofP^l\ Finally, we summarize our work in sec. 5. 



2 Recalculation of the LO splitting function 

As a first example, we will rederive the LO result for the flavour non-singlet splitting 
function, using the ML prescription. This is a rather trivial calculation that nevertheless 
displays the main improvements provided by the use of ML. Furthermore, the virtual 
graphs in the NLO calculation have the LO kinematics, so this section also serves to 
prepare the NLO calculation. We noted before that the LO example has already been 
worked out in 0] where coUinear poles were regularized by taking the initial quark off- 
shell, < 0, rather than by using dimensional regularization. This is perfectly fine at 
the LO level, but beyond LO it seems a forbidding task to keep p"^ ^ 0, and in fact the 
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underlying method of |T0[ that we are employing has been set up in such a way that 



it relies on the use of dimensional regularization, yielding final results that correspond to 
the MS scheme. It therefore seems a useful exercise to sketch the calculation of Pgg'* in 
the ML prescription if dimensional regularization is used. 



The Feynman diagrams contributing to Tgg at LO are shown in Fig. 1. For the gluon 
polarization tensor in diagram (a) we need to insert the two parts of the ML discontinuity 
in (^) with their two different 5-functions. The first part of the phase space, resulting 
from 6(1'^), can be written as 



71^-" 



driikir, (9) 



. . . 2r(l-e) 
where 

kl = ll = \e\(l-x) . (10) 
The ^(/^) contribution of graph (a) to Tgg is then given by 

= ^pp rd\e\\kr'-' cAi-xr^-±^ . (n) 

2n Jo 1- X 

Using the identity 



-i-e _ l^^i , 1 . f\n{l-x] 



{I - x)-'-^ ^ —6{1 - x) + ^-e ^ ^ +0{e'), (12) 

e [i — x)+ \ i — X J ^ 

where the 'plus'-prescription is defined in the usual way, one readily obtains 
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r(").'5('')(x) = ::^pp / d\k''\\kT''' Cf 

27r ./o 



For the ghost-like part we introduce the variable k as 

/ I — / — I /c I /t^ . 

The phase space is then given by 



(13) 



(14) 



71 



1-e 



2r(i - 6) 

-l-e 



TT 



2r 1 



-5(1 -x) 



d\k^\ / '5((1 -x)(l - k)) 

2| 







dn 







1 



(15) 



where the last line follows since the root of the delta function for k = 1 never contributes 
when we insert the second term in (^ for the gluon polarization tensor into graph (a), 
thanks to the factor 2n*l/l']_ ~ (1 — k)//? accompanying the 5{P + Z^.) in (§)• Thus, the 
ghost part contributes only at x = 1. The contribution ~ 1/k of 2n*l/l']_ gives rise to a 
1/e-pole in the final answer: 



i qq 



a 



5{l-x)VV / d\e\\k: 
2vr Jo 



2\-l-t 



Adding Eqs. (|13]) and ([16|) , we get the full contribution of graph (a) to Vq^. 

a. 



^ 1 I ^2 

d\k^\\k\^-'C, ^ 



(16) 



(17) 



'0 il-a;)+ 
An important feature of this result should be emphasized, as it will also be encountered at 
NLO: the integrand in ([T7|) is completely finite, in a distributional sense. In other words, 
using the ML prescription, we have verified the finiteness of the LO 2PI kernel q ^ qg 
in the light-cone gauge. We point out, however, that the finite 2PI kernel arises as the 
sum of the two singular pieces in Eqs. ([13D , (|T6|) . This is again a finding that will recur 
at NLO: the full discontinuity @) of the gluon propagator in the ML prescription has a 
much milder behaviour than the individual contributions to it. 



It is instructive to contrast the result in (|T^) with the one obtained for the PV pre- 
scription 0: 

/ d\e\\kX^'^CF - + 



2^ 



'1 



+ 2/o5(l - x) 
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where 



u 



rdu ^ —ln5 



(19) 



,0 u^ + S^ 

Thus, the 2PI kernel for the PV prescription has a divergent coefficient of 5(1 — a; 
resulting from the gauge denominator 1/nl and being regularized by the parameter 6. 



The calculation of the LO splitting function is completed by determining the endpoint 
contributions at x = 1, corresponding to Zg in (§) and given by the graphs in Fig. 1(b). 
They can be straightforwardly obtained^ using the UV-singular structure of the one-loop 
quark self-energy, determined for the ML prescription in [|T8]. One finds ||^: 
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1 



Up — 

27re 2 



that is. 



f{0) 



(20) 



Putting everything together, one eventually obtains 



l + x2 3^,^ 



(21) 



in agreement with . We finally note that of course the same final answer is obtained 
within the PV prescription: the singular integral Jq in (p!8|) is cancelled by the contribution 
from Zg, since we have 



7PV 



27re 



--2/o 
2 



:22) 



Thus, to summarize, the advantage of the ML prescription at the LO level mainly amounts 
to producing truly finite results for the 2PI kernels, as required for the method of JIT], [1^ . 
Furthermore, there is no need for introducing renormalization constants depending on 
additional singular quantities like Jq that represent a mix-up in the treatment of UV and 
IR singularities. 



^Alternatively one can obtain the contributions from the requirement Pqq\x)dx = [ p^ , ^. 
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3 The calculation of the flavour non— singlet splitting 
function at NLO 



At NLO, there are two different non-singlet evolution kernels, P~'(^) and P~^'^^\ governing 
the evolutions of the quark density combinations q — q and q + q — {q' + q'), respectively. 
The two kernels are given in terms of the (flavour-diagonal) quark-to-quark and quark- 
to-antiquark splitting functions by (see, for instance, Ref. fl^ ) 

p±,(i) ^ pVM ± pni) (23) 

where the last splitting function originates from a tree graph that does not comprise 
any real-gluon emission and is therefore free of any problems related to the use of the 
light-cone gauge. Thus, we do not need to recalculate P^q^^^ ■ The Feynman diagrams 
contributing to P^q''^^ are collected in Fig. ||. We have labelled the graphs according to 
the notation of |T^. We also show the graphs contributing to Zg at two loops. We will 
not calculate all of these, since this is not really required. Their role will be discussed in 
subsection 3.4. 

3.1 Virtual— cut diagrams and renormalization 

Many of the diagrams in Fig. ^ have real and virtual cuts, as has been indicated by the 
dashed lines. Let us start by discussing the contributions from the virtual cuts in graphs 
(c),(d),(e),(f),(g). It is clear that these essentially have the LO topology in the sense 
that there is always one outgoing gluon (momentum /), to be treated according to the 
ML prescription as discussed in the previous section. We recall that this means that 
there are two contributions for this gluon, one at = 0, and the other with P + = 0, 
corresponding to the gluon acting as an axial ghostQ. This immediately implies that we 
will have to calculate the loop integrals for these two situations. In addition, it is clear 
that the ML prescription also has to be used in the calculation of the loop itself, not just 

^In the next subsection we will see that for diagrams (d),(f) there are also other contributions at 
P + l'^ = 0, not just the one from the axial ghost going into the loop. However, the integration of those 
contributions proceeds in exactly the same way as outlined here. We postpone the discussion of all the 
contributions at + = for graphs (d),(f) to the next subsection. 
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Figure 2: Diagrams contributing to Tqg at NLO. 

for the treatment of the external gluon. For instance, the gauge denominator l/(r ■ n), 
where r is the loop momentum, is subject to the prescription (|I|). In short, we will need 
several two-point and three-point functions with and without gauge denominators like 
l/{n ■ r), and for both 1^ = and P + = 0. 

We point out that important qualitative differences with respect to the PV prescription 
arise here: in the PV calculation one always has = for the outgoing gluon in the 
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virtual-cut graphs, and there is no exphcit dependence on For instance, the way to 
deal with the self-energies in graphs (f),(g) in the PV prescription is identical to their 
treatment in covariant gauges: one calculates them for off-shell renormalizes them, and 
eventually takes the limit — * 0. In this way, almost all contributions of the diagrams 
will vanish since all loop integrals have to be proportional to (e < 0) on dimensional 

grounds. Only the contribution from the MS counterterm remains |[19| , |^ because this is 
the only quantity not proportional to (Z^)^*^. In contrast to this, in the ML prescription 
sets an extra mass scale. For graph (f), one therefore encounters terms ~ (^^)~'^, but also 
terms of the form ~ [al"^ + bl'^J"'^. The latter terms yield non-vanishing contributions to 
the virtual-cut result even at P = 0. This is still not the case for graph (g) since here the 
pure quark loop of course does not contain any light-cone gauge denominator and thus 
does not depend on Nevertheless, one gets a contribution from the quark loop in (g) 
for + = 0, i.e., when the gluon running into the loop is an axial ghost, corresponding 
to the second part of the ML discontinuity in (^. 

As expected, in the actual derivation of the loop integrals the property of the ML 
prescription to allow a Wick rotation is of great help. Nevertheless, some of the integrals 
are quite involved, since the ML prescription introduces explicit dependence of the loop 
integrands on the transverse components due to the identity 2{nr){n*r) = + r^. 
Furthermore, since we are interested in calculating also the contributions at x = 1, we 
need to calculate the loop integrals up to 0{e) rather than 0{1). The reason for this 
is that very often the final answer for a loop calculation with = will contain terms 
of the form (1 — (a = 1,2), to be expanded according to Eq. (|T2|). As a result, 

a further pole factor 1/e is introduced into the calculation, yielding finite contributions 
when multiplied by the 0{e) terms in the loop integrals. A similar thing happens in the 
loop part with + = 0. Here, an extra factor 1/e can be introduced when integrating 
this part over the phase space in (|T5|). The higher pole terms created in these ways will 
cancel out eventually, but not the finite parts they have generated in intermediate steps 
of the calculation. The detailed expressions for the loop integrals in the ML prescription 
are given in Appendix A. 
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For the renormalization of the loop diagrams, one needs to subtract their UV poles, 
which is achieved in the easiest way by inserting the UV-divergent one-loop structures as 



calculated for the light-cone gauge in the ML prescription in g, ^ All structures 
have also been compiled in 0. The ones we need for the non-singlet calculation are 
displayed in Fig. ^. One notices that, as expected, the structures are gauge-dependent 
and Lorentz non-covariant. Even more, the expressions for the non-Abelian quantities 
Ilj^^ and in Fig. |^ are non-polynomial in the external momenta, owing to terms ~ 
l/[r2Z]. It is an important feature of the ML prescription that these non-local terms exist, 
but decouple from physical Green's functions [Q thanks to the orthogonality of the free 
propagator with respect to the gauge vector, n^V^^{l) = (this has actually been an 
important ingredient for the proof 0] of the renormalizability of QCD in the ML light- 
cone gauge). Thus, the non-local pole parts never appear in our calculation. This is in 
contrast to the PV prescription, where one has |^ contributions from the renormalization 
constants to the calculation that explicitly depend on the external momentum fractions 
X, 1 — X. 



3.2 Real— cut graphs 



Let us now deal with the real cuts. One way of evaluating these is to integrate over 
the phase space of the two outgoing particles with momenta li,l2 (for the notation of 
the momenta, see Fig. 0), in addition to the integration over the 'observed' parton with 
momentum k. This is the strategy we have adopted for all diagrams contributing to the 
Cp and the CpTf parts of the splitting function, i.e., graphs (b),(c),(g),(h). For graphs 
(d),(f), we found it simpler to use a different method, as will be pointed out below. 

If the two outgoing partons are gluons, their phase space in the ML prescription splits 
up into four pieces, as we discussed in Sec. 1. It is possible to write down a phase space 
that deals with all four parts. We leave the details for Appendix B. 

Upon integration of the squared real-cut matrix element for a diagram, each of the four 
parts of phase space gives highly divergent results, but their sum is usually less singular. 
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Figure 3: UV-divergent one-loop structures as obtained in the light-cone gauge, using the 
ML prescription. The indices i,j (a, 6) denote quark (gluon) colours; are the generators 
of SU(3). 

This is similar to the pattern we found at LO. For instance, the phase space integration of 
graph (b) of Fig. 2 (before performing the final integration over \k'^\) is expected to give 
a finite result, since the graph is 2P1 and possesses no virtual cut. This indeed turns out 
to be the case, but the individual contributions to (b) by the four parts of phase space all 
have poles ~ l/e^ and ~ 1/e which cancel out when combined. A similar cancellation of 
higher pole terms happens for graph (h) (which of course is not finite by itself, but has 
a left-over 1/e singularity, to be cancelled by the contribution from diagram (i)). Here 
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even poles ~ occur at intermediate stages of the calculation. For those graphs that 
also have virtual cuts, the situation is in general even more complicated, as cancellations 
will occur only in the sum of the real and virtual cuts. An example for this case will be 
given in sec. 3.3. 

For graphs (d),(f), the phase space integrals become extremely complicated. This is 
due to the extra denominator l/(/i + h)"^ present in these graphs, which causes great 
complications in the axial-ghost parts of the phase space. For (f), we found it still 
possible to get the correct result via the 'phase space method', but for (d) this seemed a 
forbidding task. It turned out to be more convenient to determine the result in a different 
way: if one calculates, for instance, the gluon loop in graph (f) for an arbitrary off-shell 
momentum / going into the loop, the imaginary part of the loop will correspond to the 
real-cut contribution we are looking for. To be more precise, the strategy goes as follows: 
we calculate the loop graph for off-shell / and insert the result into the appropriate LO 
phase space. The latter can be derived as in (^, omitting, however, the there. One 
finds 



X I (fk I (fl5\x-^5{p-k-l) = —- -/ d\e\ I rfr(A;i)-% (24) 



pn J ' 2r(l 







where now 

X 

The limits for the r integration in (^) span the largest possible range for r, given by the 
conditions Z^. > 0, + /^^ > 0. The full imaginary part arising when performing the loop 
and the r integrations has to correspond to the sum over all cuts in the diagram. One 
encounters discontinuities from the following sources: 

(A) from the loop integrations. Here imaginary parts arise, for instance, if for certain 
values of r and of the Feynman parameters ti, . . . ,tk, one finds terms of the form 

/(ti,...,4,r))"', (26) 



where / is negative. Details for integrals with such properties are given in Ap- 
pendix C. The imaginary part originating in this way essentially corresponds to the 
cut through the loop itself, i.e., to the real-cut contribution we are looking for. 
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(B) from the propagator 1/iJ? + ie) via the identity 



PV( ^ ) -Z7r5(/') , (27) 



P + ie \l\ 

where PV denotes the principal value. The imaginary part ~ S{P) obviously repre- 
sents the loop contribution at = which we have determined in the last subsection. 
Therefore, we do not need to reconsider this part of the discontinuity. 

(C) from terms ~ for which a relation similar to (|27| ) holds, 

^ ^ FVi \] -iTTsign{nl)5{nl) . (28) 



[nl] nl + iesign{n*l) V^^. 

At first sight, one might think that the discontinuity ~ 6{nl) simply corresponds 
to the calculation of the gluon loop for the case when the gluon entering the loop 
is an axial ghost with P = —l'j_. However, the situation is more subtle: The terms 
~ do not only originate from the propagators of the external gluons, but also 

from splitting formulas like 

I = ±(^+ ' ) (29) 

(where r is the loop momentum), as well as from the loop integrals themselves, like 
in the case of 

d'^r 1 



(2vr)rf (r2 + ie){{l - r)2 + ie)[nr] ' 
All these terms ~ have to be treated according to the ML prescription, i.e.. 



give rise to discontinuities ~ 6{nl) ~ 6{1 — x) via (^81) . The sum of all discontinuities 
arising in this way actually has to correspond to the 'pure-axial-ghost' part of the 
graph, given by (a) the virtual-cut contribution when the gluon going into the loop 
is an axial ghost, plus (b) the real-cut contribution when both final-state gluons act 
as axial ghosts. These two parts cannot easily be separated from each other, which 
is the reason why we postponed the whole treatment of graphs (d),(f) at P = —l'j_ to 
this section. The integrals needed to obtain this part of the discontinuity are those 
already mentioned in the last subsection and collected in the right-hand column of 
Tab. 4 of Appendix A. 
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It is also worth mentioning that despite the fact that graph (f ) has a squared gluon propa- 
gator, there are cancellations coming from the algebra in the numerator; as a consequence, 
one never encounters expressions like or 1/ {t^ + ieY before taking the discontinuity. 



and (|28| ) , (|2^) are all we need. 



Clearly, when finally collecting all the imaginary parts from (A) and (C), the PV-parts 
in (B) and (C) play a role in the calculation. While 1/nl ~ 1/(1 — x) in ( pSj) only diverges 
at the endpoint at x = 1 where it is always regularized by factors like (1 — x)~'^, the 
propagator 1//^ in (|27|) in general has its singularity inside the region of the r integration: 
from (p5D one finds that > for r < 1, but < for r > 1. The principal value 
prescription]^ in ( pTj) takes care of the pole at r = 1 and leads to a cancellation of the 
positive spike for r ^ 1~ and the negative one for r 1^, resulting in a perfectly 
well-defined finite result. 



The vertex graph (d) can be treated in a similar fashion as (f). Here one calculates 
the full vertex for = 0, /c^ < 0, but arbitrary and determines the imaginary parts 
arising with respect to . This corresponds again to point (A) above, and calculational 
details are also given in Appendix C. The imaginary part from (B) is again related to the 
virtual-cut contribution at = that we already calculated in the last subsection. The 
discontinuity from (C) needs to be taken into account, and as before it corresponds to the 
full 'pure-ghost' contribution (virtual-cut and real-cut), residing at a; = 1. 

A final comment concerns graph (i). Its contribution to Tqq is essentially given by a 
convolution of two LO expressions, each corresponding to Fig. 1(a), keeping however also 
all finite terms in the upper part of the diagram, including the factor {l — x)~'^ from phase 
space (see (|n]),([g)): 

1 + ^2 ^J\yi{1-z) 



(i) 



-el 



where 



z 

^ dz 
z 



1 

® - 


■ 1 + ^2 - 


e 





/(f).(.). 



(31) 



(32) 



■^To avoid confusion, we emphasize at this point that the principal value for 1/P in ( p7| ) is well- 
defined here and not related to the principal value prescription for the light-cone denominator 1 /nl that 
we heavily criticised in the introduction. 
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Note that this is in contrast to the PV prescription where the contribution from (i) does 
not correspond to a genuine convolution in the mathematical sense. Since both of the 
convoluted functions in (|3lD contain distributions, the convolution itself will also be a 
distribution. The evaluation of (^) is most conveniently performed in Mellin-moment 
space where convolutions become simple products. Some details of the calculation and a 
few non-standard moment expressions are given in Appendix D. 



3.3 Final results 



We now combine the results of the previous subsections. The first observation is that 
for the ML prescription all 2PI graphs, and also the difference (h— i), turn out to give 
truly finite contributions to Tgg, before the final integration over is performed. This 
expectation for the light-cone gauge [^] was not really fulfilled by the PV prescription, 
where the results for the diagrams depended on integrals like Iq in ([T9|) that diverge if 
the regularization 6 of the PV prescription is sent to zero |P, 0. The finiteness of the 
kernels in the ML prescription comes about via delicate cancellations of terms sometimes 
as singular as or even when the various real-cut (gluon and axial-ghost) and 
virtual-cut contributions are added. To give just one example beyond those already 
discussed in the previous subsection, let us discuss the contributions of graph (g) to Tgg. 
From the real-cut diagram, one has up to trivial factors 



'^'^ \27iJ Jo r(l-2e) ' \3e^ 9e 3^^^ 27 J 

l + x^ / 2 10\ 4, /ln(l-x)\ 21 + x2 

+ 7 — ]+-(l + x^){— -] Inx 

l-a;+V 3e 9/ 3^ ^ V 1 - a; /, 31-x 



(33) 



The virtual-cut graph for P = — (corresponding to the gluon being an axial ghost) 
contributes before renormalization: 

r.?-(S)^*a-.^-)ppfM^(_|,_i5_^«,_p). ,3, 

The loop with = only contributes via its renormalization counterterm as explained 
earlier. This contribution exists also for the loop with = and reads on aggregate 
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for both foop parts: 

\27rJ Jo r(l-2e) 3[e(l-x)+ ^ ^ ) [ i _ a; ' + ^ 

(35) 

When adding the integrands of Eqs. all poles cancel, and as promised the 

contribution to Tgg is finite before integration over \k'^\. 

Next, we determine the contributions of the various graphs to P^q^^\ making use of 
Eq. (^. The results are displayed in Tables 1 and 2. We see that all entries in the tables 
are completely well-defined, even at x = 1, in terms of distributions, which is a property 
that we already encountered at LO. 

The sums of the various graph-by-graph contributions are also presented in Tables 1,2. 
One realizes that many more complicated structures, like the dilogarithm Li2(x), cancel 
in the sums. Considering only a; < 1 for the moment, it is the most important finding of 
this work that the entries in the columns 'Sum' in Tables 1,2 exactly reproduce the results 



found in the PV calculations for x < 1. Since the latter are in agreement with 

those obtained in the covariant-gauge OPE calculations ^2|, we conclude that the ML 
prescription has led to the correct final result. To a certain extent, this is a check on the 
prescription itself in the framework of a highly non-trivial application. Since - in contrast 
to the P V recipe - the ML prescription possesses a solid field-theoretical foundation ^ , 
our calculation has finally provided a 'clean' derivation of the NLO fiavour non-singlet 
splitting function within the CFP method, highlighting the viability of that method. 

The next subsection will address the endpoint {5{1 — x)) contributions to the NLO 
fiavour non-singlet splitting function. 
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Table 1: Final results for the Cp part of Pqq ' on a graph-by-graph basis. 
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Table 2: Final results for the CpNc and CfT/ parts of Pqq ' on a graph-by-graph basis. 
The dots in the contributions ~ 5(1 — x) in the first row denote uncalculated pieces consisting 
of C(3) and a rational number (see App. C). 
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3.4 Contributions at x = 1 and a sample two— loop calculation 

of Zq 



In the PV calculations [1^, [T2| of the two-loop splitting functions the contributions 
~ 5(1 — x) were never directly calculated, but inferred ^ from fermion number 
conservation, expressed by the requirement 

>v,{i). 



x))dx = . (36) 



We could proceed in the same way and stop the calculation here. However, the calculation 
we have performed in the ML prescription allows us to go beyond this pragmatic approach, 
since - at least for the Cp and CpTf parts - we have always picked up the finite amounts 
of 6{1 — x) contributed by the 2PI kernels. If we now performed a calculation of C,, 



(1) 

Q 5 



corresponding to the graphs in the bottom row of Fig. |^, we would have all terms ~ 6{l—x) 
in the NLO flavour non-singlet splitting function and could check whether indeed (|36|) 
is correctly reproduced. As an example, we will go this way for the CpTf part of the 
splitting function. 

Let us first establish what we need to get for ^q^\ The coefficient of 5(1 — x) in the 
NLO splitting function was determined in via (p6| ) to be 



(37) 



C| Q - 3C(2) + 6C(3)) +C^Tf (-1 - ^C(2)) +C^Nc + yC(2) - 3C(3)) 

while in our calculation we have according to Eq. (^) and Tables 1,2: 

e^'^ + Cl(^l- 3C(2) + 6C(3)^ + CpTf (-^C(2)) + CpNc C(2) + • • , (38) 

where the dots indicate that we have not entirely calculated the finite amount of 5(1 — x) 
in the CpNc part of Pqq^\ even though we were able to determine its contribution ~ ^(2) 
(see App. C). Comparing Eqs. (|37|) and (|38D, we get a prediction for the and CpTf 
parts of ^g^^ in the light-cone gauge with ML prescription: 

e^') = -1^1- \cFTf + CpNc (...), (39) 

where all we can say about the CpNc part is that it does not contain any terms ~ C(2)- 
It is quite remarkable that no C(2), C(3) terms are left over in the C'jp and CpTf parts of 
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To directly calculate the CpTj part of ^q^\ we only have to evaluate the first diagram 
in the bottom row of Fig. ^ What we need to extract is the two-loop renormalization 
constant for that diagram, when the light-cone gauge with the ML prescription is used. 
The calculation is relatively easy since the inner quark loop has obviously no light-cone 
gauge propagator and can in fact be calculated exactly: 

^ ^ a, 8r2(2-e)r(e) / 47r \' r 2 1 

n..(r) = -^Tj^^^^^i^ [-Z^.) - V^] • (40) 

This self-energy can then be renormalized with the help of the counterterm in Fig. |^. The 
renormalized loop is then inserted into the outer loop. Here it is very convenient that II^i, 
is transverse, 



Q/3/ 



(41) 



In other words, the whole calculation is not very different from a simple one-loop calcu- 
lation of the quark self-energy, the only exception being that we now need loop integrals 



with the extra factor f— r 



2\-« 



present in 



If we embed the whole graph into the 



Dirac trace as shown in Fig. |l|(b), it turns out that we only need few integrals of this 
kind; they are collected in Appendix E. Since we have renormalized the inner loop, the 
left-over divergence after loop-integration determines the two-loop counterterm and thus 
the contribution to We find in the MS scheme: 



Comparing to Eq. 
it to be in dSQl): 



this implies that the CpTj part of ^^^^ is exactly what we expected 



(42) 



--CpTf . 
6 ^ 



(43) 



This result clearly demonstrates the consistency of the whole approach: our example shows 
that the light-cone gauge method of [0, |l^ is also able to determine the contributions 
~ 5{1 — x) to the splitting functions by explicit calculation. It would be interesting 
in this context to calculate also the other contributions to ^^^^ (and the missing part 
~ CfNc5{1—x) in our Tab. 2); important steps in this direction have been taken in [24, Q 
by examining the second and the third diagram in the bottom row of Fig. |], which yield 
the C|n part of Indeed it turns out that the results of ||2^, ^ can be exploited to 
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reproduce the term — 5C|./8 in our prediction (pOf ) for ^q^\ which can be regarded as a 
further confirmation of our results. 

We have to admit, however, at this point that the abihty to obtain the correct endpoint 
contributions is not restricted to the ML prescription: this is also possible for the PV 
prescription. Here, the coefficient of 5(1 — x) in the CpTf part of P^q^^^ reads 

e^Jv^^^- - C,t/-^Io . (44) 

Here the second term originates from the entry ' — 10/9' in the last column of Tab. 2, when 
we omit the 'plus'-prescription there and reintroduce it using the PV identity 

^ 1,6(1 - x) + —^-^ , (45) 



1 — X (1 — X 

where Jq is as defined in ([l9|) . Furthermore, ^l^pv^^''^ denotes the CpTj part of ^q^^ when 



the PV prescription is used. The explicit calculation gives 



2 



Z^ll' = 1 + ( ^ ) C^T. 



that is 



(46) 



Cv'" = C,T, (-1 - ic(2) + f /„) . (47) 

It is interesting to see how upon combining Eqs. ( |iiD and ( |47| ) the Jq terms drop out, 
and the CpTf part of the endpoint contributions comes out correctly as in (|37D also for 
the PV prescription. We note, however, that again this happens at the expense of having 
renormalization constants depending on singular quantities like Jq that represent a mix-up 
in the treatment of UV and IR singularities. 
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4 The calculation of the N'^ part of the singlet 
splitting function Pgg at NLO 



Let us now turn to the calculation of Pgg' . We restrict ourselves to its N"^ part, since 
the contributions ~ CpTf, NcTf are essentially trivial as far as the treatment of the 
LCA is concerned: The CpTj part comprises no gluon emission at all, and all diagrams 
contributing to the NcTf part contain a quark loop and the emission of at most one gluon. 
Such diagrams with one-gluon emission have the LO kinematics and will not reveal any 
new features as compared to what we have already discussed. In contrast to this, the 
part of Pgl^ requires the renormalization of the non-Abelian part of the three-gluon 
vertex and therefore really provides a further challenge for the ML prescription. 

The diagrams contributing to the N^. part of Tgg at NLO are shown in Fig. 4. We do 
not show here the graphs contributing to Zg at two loops, since we will not attempt to 
calculate them. 




Figure 4: Diagrams contributing to the A^^ part of F^^ at NLO. 
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The calculation of the various real-cut and virtual-cut diagrams proceeds in exactly 
the same way as before. For the renormalization of the virtual-cut contributions in the 
triangle graph (d) and the 'swordfish' ones (sl),(s2), we need the UV counter-term for 
the three-gluon vertex in the light-cone gauge with ML prescription. Here we can rely 
on the result presented in (see also 0); the part of it that is relevant for our 
calculation is recalled in Fig. 5. 



Concerning the real cuts, we mention that for graphs (h),(b),(j),(k) we use the expres- 
sion in (|B.5|) for the phase space. As for the CpNc part of P^q^^\ we found it easier to 



determine the contributions of the real cuts of the remaining diagrams via the extraction 



Pi,aif 3 III 





P2, 0-2 



perm. 



^/<u.2M = Smsm (P2 - P3)„i + 5„3/« (P3 - + fi',.1,.2 (Pi - Pa),, 
+5wf2<3 (Pi -P2) -n 



sum r^.i/.2/.3(PbP2,P3) = T^iW2t^3{Pl,P2,P3) + S^,,f,,^,,{pi,p2,P3) 



_ iga, Nc faia2C 
Air e J 



11 A — on -\- 

3 ^/ii/i2M3 IJ.1IJ.2IJ.3 ^ 



Figure 5: UV-divergent structures of the non-Abelian part of the three-gluon vertex as 
obtained in the light-cone gauge, using the ML prescription. pi,p2,P3 denote the momenta 
of the external gluons, 01,02,03 are the associated colour indices (/"i"2a3 being the structure 
constants of SU(3)), and yUi,/i2,/i3 are Lorentz indices. The dots indicate structures (some of 
them non-local) which do not contribute to our calculation thanks to the orthogonality of the 
free propagator to the gauge vector n. 
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of the imaginary parts of the associated virtual graphs. 



We have verified that again for the ML prescription all 2PI graphs give truly finite 
contributions to Tgg, before the final integration over \k'^\ is performed. This also applies 
to the endpoint x = 1, where the result for each graph is again always well-defined in 
terms of distributions and, as before, also has a coefficient of 5(1 — x) that contains no 
1/e poles. Table 3 presents the contributions of the various diagrams to P^l\ Here we 
have defined the functions 

(1 — X + x^)^ 



x{l — a;)+ 



U.) . ,48) 

/±-\-x (12' /I P' \ 1 

— ln( -] =-2Li2(-x)-21nxln(l + x) + -ln2x-C(2) . 
^ z \ z / 2 

We mention in passing that graph (j) and the 'swordfish' diagram (si) give vanishing 
contributions to Pgl^ if the PV prescription is used, but are non-vanishing for the ML 
prescription, where finite contributions arise from their ghost parts. 

As for the case of P^^^^\ the full result for the part of given by the column 
'Sum', is (at x < 1) in agreement with the PV result of [1^, which in turn coincides 



with the OPE0 one [|l^. Thus, the CFP method with ML prescription has also led to 
the correct final answer in this case, which clearly constitutes a further non-trivial and 
complementary check. As can be seen from Tab. 3, we have not determined the finite 
amounts of contributions ~ 6{l—x) for the graphs since, like in the case of the CpNc part 
of P^q^^\ these are quite hard to extract in some cases. The endpoint contributions to 
Pgl^ can then only be derived from the energy-momentum conservation condition [12 . 



We emphasize however that, just as for P^q^^\ there is no principal problem concerning 
the calculation of the endpoint contributions: had we calculated the full 5(1 — x)-terms 
in Tab. 3 and the two-loop quantity ^^g \ all endpoint contributions would be at our dis- 
posal, and it would no longer be necessary to invoke the energy-momentum conservation 
condition; in fact, this could serve as a further check of the calculation. 



See also our discussion in the introduction concerning the OPE calculations |l5m of 
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5 Conclusions 



We have performed a new evaluation of the NLO flavour non-singlet splitting function 
and of the Nq part of the NLO gluon-to-gluon splitting function within the light-cone 



gauge method of |Tl|, . The new feature of our calculation is the use of the Mandelstam- 
Leibbrandt prescription for dealing with the spurious poles generated by the gauge denom- 
inator in the gluon propagator. In contrast to the principal value prescription employed 
in previous calculations |10|, , the ML prescription has a solid field-theoretical foun- 
dation and will therefore provide a 'cleaner' derivation of the result. As expected, the 
final results come out correctly, i.e., are in agreement with the ones in [l^, |2^ . 
This finding is both a corroboration of the usefulness of the general method of |]TT], ^ to 
calculate splitting functions, and a useful check on the ML prescription itself in a highly 
non-trivial application. 

We have also discussed the 5{1 — x) contributions to the NLO flavour non-singlet 
splitting function, performing an explicit sample calculation of a two-loop contribution 
to the renormalization constant Zq in the ML light-cone gauge. It turns out that one 
indeed obtains the right amount of contributions at a; = 1 as required by fermion number 
conservation. 

We conclude by conceding that the ML prescription is in general much more compli- 
cated to handle than the simpler, but less well-founded, PV prescription. With regard to 
future applications at, for instance, three-loop order, this creates a certain dilemma: the 
ML prescription might be too complicated to be used in that case, while on the other hand 
the ill-understood success of the PV prescription at the two-loop level is not a warranty 
that it will also produce correct results beyond. 
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Appendix A: Virtual integrals 



Here we list some loop integrals needed for the calculation. We do not need to recall any of 
the covariant integrals, which are standard, but will only present those with a light-cone 
gauge denominator, to be treated according to the ML prescription (|l|). 

We begin by performing a sample calculation of the integral 

d'^r 1 



I{n,q) = 



{27tY (r^ + ie){{q — r)^ + ie)[nr] 
d'^r n*r 
(2tt)^ (r^ + is){{q — r)^ + ie){nrn*r + ie) 



(A.l) 



We recall the definitions 



n=^{l,0,...,0,-l) , n* = —{1,0,..., 0,1) = —p 



(A.2) 



where p = P(l, 0, . . . , 0, 1) is the momentum of the incoming quark, see Fig. |^. Introducing 
Feynman parameters, one has 

d'^r 



4P /"^ Z"^"* 
I(n, q) = — dt ds 
pn Jo Jo 



To - 



(27r) 



+ sr\_ — 2(g ■ r)t + qH + ie 



(A.3) 



After performing a Wick rotation and straightforward integrations over r one arrives at 
iV{l + e) 2n*q ( An 



I{n,q) 



q2 _|_ V_g 



dtr'(l-t)-^-' [ ds(l + st, ^ . 



q +qi + le 



For example, for the case q = k one finds 

zFfl + e) / 47r \^ 1 



I{n, k) 



C(2) - Li2 



+ 2eC(3) 



(A.4) 



(A.5) 



where we have kept those terms that contribute to the final answer. In ( |A.5D , ({n) is 
Riemann's C^function and Li2(a;) denotes the dilogarithm, defined by 

_ f' ln(l - zy) 

Jo y 



-dy . 



(A.6) 



The result in ([A.5|) coincides with the one in for e = 0. Note that the ML prescription 
arising for l/[nA;] is actually immaterial here since nk = xpn never vanishes. 



28 



Setting, on the other hand, q = I one gets for 

XI + e) / 47r 

167r2 \-\k^\(l-x)J nl 



/(„,0 = £il±£lf^J^^)'lJ:,, (A.7) 



and for = — 



Note that the real part of (|A.7| ) has to be takenQ. Table 4 contains all the required 
integrals with an ML light-cone gauge denominator. The integrals in the first column are 
for /2 = 0; they depend on 

nk _ nl /AX 

X = — and X = — = 1 — a; . (A. 9) 

pn pn 

Recall that terms ~ (1 — a;)~^~"^ will lead to further poles, as was shown by (p!2D. The 
integrals in the second column of Tab. 4 are for the axial ghost case, P = —l^, and 
eventually need to be integrated further over the variable k defined in Eqs. (|T3p,(|np. The 
K integration produces further poles. We have 

K = — ^7 and K = 1 — K . (A. 10) 

F I 

We note that the last integral was much easier obtained by performing the k integration 
before the ones over the Feynman parameters. Therefore we only present the final, k- 
integrated, result in this case. As can be seen, the integral was accompanied by two 
different powers of k. 



^Here one obviously has to discard the overall factor i. 



29 







f + f, 







pn 



r'^{k—r)'^[nr\ 



C(2) - Li2 (x) + 2eC(3) 



C(2) - Lis («) + 2eC(3) 



(r^— 2nfc n*r) pn 
r'^{k—r)'^[nr\ 



i + 2(l-lnx) + lLi2(x) 



-iC(2) + e(4 + C(2)-2C(3)) 



i + 2(l-lnAi) + iLi2(/€) 



4C(2) + e(4 + C(2)-2C(3)) 



pn 



^i;-i-(l-e2C(2) + 2e3C(3)) 



-e 1 



2e2 



(1-^) (l + e'C(2)) 



(r^—2nl n*r) pn 
r2(«-r)2[nr] 



-2A;^5 



i + 2 + 6(4- C(2)) 



pn 

'(p^^r)^(fc^^rp"[nrj 



iif - iln2a;-Li2(x) 



-e5 + 2eC(2)lnx 



(i + 2 + 4e + eC(2)) - 1) 
+ (l-i) ln^ + 6(3-C(2)) 



pn 



(p— r)2(i— r)2 [nr] 



does not occur 



pn 



(k+r)'^(l~r)^ [nr] 



[X 



x-)(^-C(2) + 2eC(3)) 
-3C(2)x-^ 



i + 2C(2) + 2eC(3) 



pn 

r'^ {p—r)'^ (k—r)'^ [nr] 



^ + ^ - 2Li2(x 
iln2x-2eC(3) 



^ (1 - '^-^) 

-/t-^ (Li2(/€) + C(2)) 

-2 (Li2(/i)-C(2) + eC(3)) 



pn 



r2(fe+r)2(i— r)2[nr] 



A:2 



3 + in£ _ Li2(a;) 



2e2 



1 



-iln2a;-|C(2)-e(5 + 5C(3)) 
■|elnx(ln^x — 2Li2(x) — Inxlnx) 



i In K + Li2(/t) — I In^ K 
K-(i + 3C(2) + 6eC(3)) 



pn 



r^"(p^^r)^(I— r)^"[nrj 



does not occur 



VI 



2e 



' ^ - 3C(3) 



(after integration dnn 

f[7 + 3 + C(2); 
(after integration jj^ dnn^' 



Table 4: Two- and three-point integrals with a light-cone gauge denominator for 
the ML prescription, calculated up to C(e). We have dropped the ubiquitous factor 
i/167r2 (47r/|P|)'r(l - e)/r(l - 2e). x and k have been defined in Eqs. ( |OD and (|CTD , 
respectively. 
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Appendix B: Three— particle phase space 



As we discussed in Sec. 1, the phase space for two gluons (plus one 'observed' parton) 
will split up into four pieces for the ML prescription: 



PS2 = X j d%d% (5 (^1 - X - Si^l + ll±Ml) , 

PS, = xj d\d% 5[l-x- ^{llWl + 111) , 

]di .di ^ nli + nl2\ ^,j2 , ,2 ^.n2 , ,2 



PSi = X I d%d% 6{l-x I + K±)^i^2 + ^2,±) , (B.l) 



where li,l2 are the gluon momenta. The 5 functions in ( p.l|) determine whether one (or 
both) of the gluons acts as an axial ghost. 

As we know from the discontinuity in (|^), the tensorial structures of the non-ghost 
part and the ghost part are different. However, we can rewrite (H) as 



971*] 

disc [P^"(/)] = 27re(/o)-^ ( 5(/2) - 6{P + ll) 



(B.2) 



This is possible because of 2n*l/l\ = 1/nl for = and {nl){n*l) = for l"^ + !]_ = 0. 
In this way, it is always possible to calculate just one combined matrix element, using 
the tensorial structure in square brackets, and integrate it over a phase space subject to 
simply the difference — 5{f' + For our two-gluon case, this means that we have 
to consider only the combination 

PSi - PS2 - PS3 + PS4 . (B.3) 

We now introduce the Sudakov parametrizations 

pn 

^2 = z(l-y)pf^+^^n^ + li^^, (B.4) 
pn 

where (^f_|_)^ = —^i±- The first 6 functions in ([B.l| ) imply y = x/z. If one wants to 
integrate over an arbitrary function / of scalar products of the momenta, one writes the 
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four parts of phase space in the following way: 



nli+nl2\ ,fll± -{p-kf hp 



xld%d%5{l X J/V|F|' |P| 'IPI'l^i'-^' 



X, z 



1 Tl^-' 



7r2 



4r(l-e)r(i-e) io io io 



de sin"2. ^ 



/2 



a^ ''f\—vw{l—vx) 



XV 



X 



\k'^\vx 1 — vx 



- {rl + r\ — 2rir2 cos 6*) , 



'■2,. 



2|A;2|xf 



X, 1 — f X 



X 



^/(l!!^i,™,i(.| + .J-2.r,r4Cos«),^,x,l 



f V X X 



2\k^\x 



\ ' 5 ' 5 

V X X 

1 



r? + Tg — 2r^rQ cos 9, 



2x 



where 



+^ / ■ rfpp--p-</(^. -. r? + r| - 2r,r,cos0. 1. 1 

X = 1 — X , V = 1 — V , w = 1 — w , p = 1 — p 
ri = , r2 



(B.5) 



vwvx"^ 



X 



r5 



w , r4 



vwx 



wx 

X 

wp 




re 



X 



'vw 



vw . 



(B.6) 



The meaning of the arguments of the function / in the various parts of the phase space 
can be seen from the first line of ( |B.5| ). We note that the last part of the phase space 



in (p.5| ) is proportional to 5(1 — x); it corresponds to the contribution where both gluons 
are axial ghosts. 



We do not list the results for the numerous different phase space integrals one encoun- 
ters. We just mention that one frequently needs the integral 



de 



sin-2^ e 



r1 + Ti — 2rarb cos 9 



2 '2 



1 



+ ^2^1 1,1 + e,l-e;^ e(r2-r2 



(B.7) 
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where 2F1 (a, b, c; z) denotes the hypergeometric function [p3] . To expand in e, the relations 



2Fi(l,l + e,l-e;z) = (I - z)-'-'' 2Fii-e, -2e,l - e; z) , 



(B.8) 



2Fi(ae, be, I - ce; z) = 1 + abe"^ (Li2{z) + ecLi3(z) + e{a + b + 0)81,2(2;)) + 0{e' 



are useful, where 



LUz) 



' Inylnjl-zy) ^ , 1 f' \n\l-zy) ^ 
dy , Si,2 2:) = - / dy 

y ^Jo y 



(B.9) 
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Appendix C: Imaginary parts of loop integrals 



To see how to extract the imaginary part of a loop integral, let us go back to our example 
in (A.4) for the case q = I there. The integration over the Feynman parameter s in (A. 4) is 



trivial and can be done immediately. Rather than performing straightaway the integration 
over t to get the general result of the integral for arbitrary l^, it is more convenient to 
include the r integration of in the calculation and carry it out first: 



(C.l) 







X 



where we have used the definitions in (^). To result in an imaginary part^, the limits 
for the T integration have to be chosen in such a way that those terms in (|C.1|) , that are 
raised to the power — e, become negative, i.e., < r < 1/(1 — tx) for the first term in 
square brackets in ( p.l|) and < r < 1 for the second. The r integrations become trivial 
then and lead to simple beta-functions. Afterwards, the t integration can be performed; 
the result is given in Tab. 5 where we also list other integrals that we encountered. As 
can be seen from Tab. 5, we also needed some integrals with an extra factor r or (1 — r) 
in the numerator. We do not consider the covariant integrals in Tab. 5 since in their case 
the extraction of the imaginary part is rather straightforward. 

As we discussed in Sec. 2.2 (see Eq. (p7|)), we have terms ~ 1/P ~ 1/(1 — r) in the 
calculation, resulting from the propagator of the gluon running into the loop, and to be 
treated according to the principal value prescription. Therefore, we will also need integrals 
like ( |C.1|) with an extra factor 1/(1 — r) in the integrand. Such integrals are in general 



much more difficult to calculate. A typical integral needed is 



PV 



L "'0 



1 - T 



ml^iil. (C.2, 



The integrals with an extra denominator 1 — r are also collected in Tab. 5. 

We finally emphasize that Table 5 does not contain all terms ~ e, but is in general 

correct only to 0{1). The only 0{e) terms that are fully accounted for in Tab. 5 are 

^We obviously do not take into account the overall factor i here. 
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those ~ C(2)- We have not consistently determined the other contributions ~ e, because 
this is a very hard task. Therefore, since terms hke (1 — will lead to further pole 

factors ~ (5(1 — x)/e via Eq. ([l^), we will not be able to calculate the finite amount of 
5(1 — x) in the final result for the CpNc part of the two-loop splitting function, except 
for the contributions ~ C(2)5(l — x). However, the expressions in Tab. 5 are sufficient for 
checking graph by graph the cancellation of all pole terms proportional to 5(1 — x), i.e., 
for proving the finiteness of the 2PI kernels at x = 1 in the ML prescription. 
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-XV 



pn 



r'^(l—r)^lnr] 



J drr-^r°(l-r)'3 

(l + 2e)i^ + eC(2) 
(a = 0,p = 0) 



+i(Li2(x)-C(2)) 



{r^—2nl n*r) pn 
r2(«-r)2[nr] 



does not occur 



— 1 — i + In X 



lnx(l - ^) - 2x- 



pn 



(p— r)^(fc— r)^ [nr] 



—X ^ ^Inx 
(a = 0,/5 = 0) 



^ + iln^x 



-2Li2(5) + e£C(2) 



(p— r)^(i— r)^ [nr] 



-i (« = 0,/? = 0) 



does not occur 



pn 



(fc+r)^ (Z— r)2 [nr] 



i + 2 + In X 

(a = 0,/3 = 0) 



(l-x-)(^-C(2))-^ 
+C(2) - 2Li2(x) - iln^x 



pn 



r^(p— r)-^(fc— r)-^ [nr] 



-i + f lnx + eC(2) 
(a = 0,/5 = 0) 



1 I 2\nx 
.2 "T ^ 



-4Li2(l-x)-C(2) 



pn 

r2(Arrr)2(I3;^)2p^ 



1 



_-(l-5-^)(^-C(2)) 
— Li2(x) — In X In a; 
(a = 0,P = 0) 
i + 2-21nx(l-i) + lnx 
(a = l,p = 0) 



1 



3™ 



5^-1^ -C(2)) 



+2(Li2(x) - C(2)) 



pn 

r2(p— r)2(Z— r)2[nr] 



a; 



-(l-x-)(^-C(2)) + ii^ 
+2Li2(x) + iln2x-3C(2) 

(a = 0,p = 0) 
^[-f -2 + 3x + 2xln5; 
(l-3x)lnxl {a = 0,(3=1) 



does not occur 



Table 5: Imaginary parts 'XV of loop integrals for the ML prescription, after integration 
over the variable r of the LO phase space (|24l). As before we have defined x = 1 — x. The 
integrals are in general only correct to 0{1); see text. Note that terms ~ (1 — x)"" must 
not be expanded in e, as further pole terms can arise via Eq. (P^. We have dropped the 
ubiquitous factor l/16n^ (47r/|A;2|)' r(l - e)/r(l - 2e). 
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Appendix D: Mellin moments 



The Mellin-moments of a function f{x) are defined by 

f dxx--'f{x). 

Jo 

As a result, the moments of a convolution f ® g (see Eq. 
moments of / and g: 

U®9T = r g^. 



(D.l) 



2|)) become the product of the 



(D.2) 



The moments of (^) are easily obtained using the formulae in the appendix of |^ . To 
invert the moments of the product f^g"' back to x-space, one needs some further moment 



expressions. Everything can be derived from the relations in ||3^, and from 



dx x"-Mn^ri -x) 



dx X 



n—l 



-x) 

1 — X 

dx x""-^Li2(a;) 

dx x"^'^ In X ln(l — x) 

^ , „_i Inxlnfl - x) 

dx X 



'0 

where 



1 — X 



-(Sl{n) + S2{n] 
n V 



ll{sl{n) + S^{n)) - ^Sl{n) - S,{n)S2{n) - \s,{n) 
-C(2)-i^5i(n), 



n 



-fo(n)-C(2))+^^iH 



(D.3) 



S,{n) - C(2) ) ( i - S,{n)] + \Sr{n) - S,{n) + C(3) 



n 



Skin) 



n ^ 



(D.4) 
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Appendix E: Two— loop integrals 



For the calculation of the CpTj part of the two-loop quark self-energy we need some 
integrals with an extra non-integer power of — in the integrand, where r is the loop 
momentum. Making use of the identities 

1 , x"-^ 

a / dx- 



a^h Jo [ax + 6(1 - ' 

a{a + l) dx dy- , (E.l) 



a'^bc ■ ''Jo ""Jo [ax + 61/ + c(l - X - ' 

one obtains rather easily: 

d'r {-r'r ^ ^ 6r(26) 1(1 - e)r(l - 26) 

(27r)<i (p - r)2 167r2^^^^^ r(l + e) 1(3 - 3e) '^^^ 



(27r)'^r2(p-r)2 167r2' ^ v ^/ r(l + e) 1(2 - 3e) 

(-^^)"^ = fV)-^fl + eV(2)) ^ (E4) 

J {27cY r\p - ry[nr] ^ ^ ^ ^ ^ J {271^ r^{p ~ ry[nr] ' ^ ' ' 

where the integral on the right-hand-side of ( |E.4| ) was determined in App. A and is 



actually finite. For the ML prescription we therefore do not need the integral on the 
left-hand-side of ( [E.4|) ; however, we will see below that the integral is divergent for the 
PV prescription. Also note that the integral in ( [E.2| ) vanishes if the factor (— r^) is not 
present. 

Finally, for the PV prescription one obtains for the integral in ( |E.4| ): 
while 

/ ftX 2f \2r 1 = T7f^(4^)^ i'P'r — [^0 + <(2) - eh] + 0{e) . (E.6) 
J (/vr)'^ r^(p — rj"^[nrj lovr^ epn 

Here 

^ i ^''^ ■ (E-'' 
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